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Abstract 

For lattice systems of statistical mechanics satisfying a Lee- Yang property (i.e., for which 
the Lee- Yang circle theorem holds) , we present a simple proof of analyticity of (connected) 
correlations as functions of an external magnetic field h, for Re h ^ 0. A survey of models 
known to have the Lee- Yang property is given. We conclude by describing various applica- 
tions of the aforementioned analyticity in h. 



1 Introduction 

Approximately sixty years ago, while studying phase transitions in certain types of monatomic 
gases, Lee and Yang were led to investigate the location of zeroes of the (grand) partition function 
of such systems as a function of external parameters, in particular the chemical potential, as the 
thermodynamic limit is approached [15]. At first sight, obtaining information on the location 
of these roots appears to be a formidable task. Yet, for a classical lattice gas with a variable 
chemical potential equivalent to the Ising model in an external magnetic field, they were able to 
show [16j that the distribution of zeroes exhibits some astonishing regularity: all the roots lie on 
the imaginary axis (i.e., on the unit circle in the complex activity plane; see Theorem [H below). 
Their celebrated result was subsequently extended to plenty of further models. 

In this note, we offer a short review of the present "status quo" regarding this so-called Lee- 
Yang theorem (Section 2), and then sketch some applications, using some novel arguments. Our 
main results can be found in Sections [3] and [4] Here is a brief summary. For simplicity, we consider 
models on the lattice With each site, x, of the lattice we associate a random variable (a 
"spin"), a x , taking values in a measure space O C R-^, for some N € N. The a-priori distribution 
of this random variable is specified by a probability measure, /j,q, on £1. To an arbitrary finite 
sublattice A CC Z d , there corresponds a space, = {<ta := {o~ x }xeA ■ ffj: £ 11, Vs £ A}, of 
spin configurations in A. Interactions between the spins are described by a potential, <£, which 
associates with every ICC^a continuous function & (X) € C(£l x ) representing the interaction 
energy of all the spins in X. Given a finite set A CC 7L d , the Hamilton function or Hamiltonian 
(with free boundary conditions) of the system confined to A is defined by 



Ht=Y,®( x )> (1-1) 



XcA 
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and the corresponding partition function at inverse temperature /3 > by 

Z ftA ($) = / e^ H tM TT d^ (a x ). (1.2) 

By £> we denote the "large" Banach space of interactions consisting of all translation-invariant 
potentials <1> satisfying |||5>||| := ^xgol^f 1 ' II^H^OIloo < °°- It is a classical result (see for 
example [26], Theorem II. 2.1, or [22], Theorem 2.4.1) that, for $ £ B and arbitrary j3 > 0, the 
free energy density 

/(/?,<&) --/T 1 lim JAI^log^AC*) (1-3) 

exists and is finite, where the thermodynamic limit is understood in the sense of van Hove. 

In what follows, the potential <3? always includes a contribution due to an external magnetic 
field, h. For example, if N = 1 the Hamiltonian is of the form 

Ha((?a) = H a ' (<7a) - h^2a x , where i?*' (o- A ) = Hf(a A )\ h=Q , a x G R. 

For certain choices of (<&,/io) (see Section [2] for an overview), the partition function ^ )A ($) is 
known to be non-zero in the regions M± = {h € C : ±Re/i > 0}. This is the famous Lee- Yang 
theorem. 

The main result of this note can be formulated as follows. Assuming that the pair (<J>,/zq) 
and the boundary conditions imposed at dA are such that the Lee- Yang theorem holds (and 
under suitable assumptions on the decay of the measure jiq, at infinity), 

the connected correlation functions (a xi ; . . . ; cr Xn )\^^h are analytic in h and ^ 
have a unique thermodynamic limit analytic in h in the Lee- Yang regions ELj-, 

where arbitrary sites in 7L d and j3 > 0; (see (|3.2p for the definition of connected 

correlations). The proof of this result is given in Section [3j see, in particular, Proposition [71 At 
first, A is chosen to be a rectangle and periodic boundary conditions are imposed at dA. But 
the result can be seen to hold for arbitrary boundary conditions for which the Lee- Yang theorem 
is valid. Earlier results of this kind can be found in |11[ [12], where the example of the Ising 
ferromagnet is treated. Our methods enable us to extend (jl.4p to systems of multi-component 
spins (N > 2) satisfying a Lee- Yang theorem, as discussed at the end of Section [3j Such systems 
include the rotor and the classical Heisenberg model, with suitable ferromagnetic conditions 
imposed. Similar results (for "Duhamel correlation functions") can also be proven for certain 
quantum- mechanical lattice spin systems. 

In Section 21 we discuss some applications of (jl.4p to classical spin systems and scalar Eu- 
clidean field theories. Assuming that the magnetic field h is different from 0, uniqueness of the 
thermodynamic limit of correlations for a large class of boundary conditions is proven, proper- 
ties of the magnetization and of the correlation length are reviewed and some bounds on critical 
exponents for the magnetization and the correlation length as functions of the magnetic field are 
recalled. 



2 The Lee- Yang theorem - a tour d'horizon 

In this section, we present an overview of classical (and quantum) lattice systems for which a 
Lee- Yang theorem is known to hold. For the sake of clarity, the models of interest are divided 
into four groups. 
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Ising-type models 

These are one-component models (N = 1) with a Hamiltonian H^, A CC Z, d , given by 

H\(a\) = - ^2 JxyVxVy — ^2h x a x , (2.1) 

{x,y}cA zGA 

where the first sum is over all pairs in A ,J XX = 0, for all x, and J xy = J yx , for all x, y. Among 
the best results concerning this class of models is one established by Newman [19j, which is 
summarized in the following theorem. 

Theorem 1. If the pair interaction in (|2.ip is ferromagnetic, i.e., the couplings J xy satisfy 

Jxy — Jyx ^ 0> Vx, y, (2.2) 

and llq is an arbitrary signed (i.e., real-valued) measure on R that is even or odd, has the property 
that f K e bcr d\fjLo(cr)\ < oo, for all b>0, and satisfies the condition 

fa(h):= / e ha dfi (a) ^ 0, V/t G M+ := {z G C : Re z > 0}, (2.3) 
Jr 

then, for arbitrary {3 > 0, the partition function corresponding to H\ in (|2,ip satisfies 

Zp,A[{h x } xeA ) = f e -^( CT A) JJ dAio((Tx ) o if /i x G M+, Vx G A. (2.4) 

The condition (|2.3p is quite natural in that it requires (|2.4p to hold for non-interacting spins, 
i.e., J xy = 0, for all x,y. (Clearly, (|2.4p trivially implies (|2.3p .) Moreover, Theorem [1] continues 
to hold if, at different sites of the lattice, different a-priori distributions are chosen, provided 
each Hoxt x G A, satisfies the conditions on no formulated above. 

In the (physically most relevant) case of positive measures, Lieb and Sokal |18| have obtained 
rather deep insights regarding the Lee- Yang property (|2.4p . In order to summarize salient features 
of their findings, we introduce some further notation. For any a > and n G N, let ^4™ + be 
the (Frechet) space of entire functions on C n satisfying ||/||& := sup zgC n e~ b ^ 2 \f{z)\ < oo, for 
all b > a. Given an open set O C C n , we define V(0) to be the class of polynomials defined 
on C" that do not vanish in O. We denote by 7 3 a +(0) its closure in -4™ + . Given an arbitrary 
distribution /x G 5'(R n ) (the space of tempered distributions on R n ), we say that \i G T n if, in 
addition, e a ' x ' G 5'(R n ), for all a > 0. A finite, positive measure fx on R" is henceforth 

called a Lee-Yang measure if fi G T n and if its Laplace transform, /t, satisfies fi G 'Po+(H+)) 
where EFJ. = M + x • • • x M+ C C n and M + = {z G C : Re z > 0}, as before. 

Theorem 2. /Tffl If p is a Lee-Yang measure on W 1 , B G P +(H!jl) for some a > 0, B is 
non-negative on supp(^) C R n and strictly positive on a set of non-zero \i-measure, then B\x is 
also a Lee- Yang measure on R n . 

This generalizes Theorem [1] (for positive hq). Indeed, first note that the product measure 
dn(a\) = Y\ xe \ dfio(a x ) with liq G T 1 satisfying ([2.3|) is a Lee- Yang measure on Rl A '; see |18| 
Corollary 3.3]. Let B be the Boltzmann factor corresponding to the Ising pair interaction in 
(|2.ip . i.e., B{as) = exp [ J2 X y( z A JxyO~ x o- y ] , with J xy G R. It is not hard to see that B G «4mjm + , 
where || J|| refers to the matrix norm of J = (J X y) x ,yeA when R' A ' is equipped with the Euclidean 
norm. One then finds |18[ Proposition 2.7] that B G V\\ (HT^ ) if and only if J xy > 0, 
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for all x,y G A. Under this condition it follows from Theorem [2] that the partition function 
Z/3,A{{h x }xeA) m Theorem [1] satisfies 

Zp = i A ({h x } xeA ) = Sii[{hx}xeA) G V 0+ [W l £), (2.5) 

which, by virtue of Hurwitz' Theorem [TJ p. 178], implies that (j2.4|) holds. Furthermore, (|2.5D 
says that 1?, the Boltzmann factor pertaining to the ferromagnetic Ising pair interaction, is a 
"multiplier" for Lee- Yang measures. 

The class of admissible measures fj, in Theorem [2] may be further enlarged by specifying 
some falloff at oo, see |18^ Definition 3.1], which amounts to requiring in Theorem [1] that 
f^e brT d\no(°~)\ < °°7 f° r some b > only. In this case, f|2.4() only holds for sufficiently small 
(3 > 0, depending on the choice of {J xy }', see |191 Remark 1.1] and |18| Corollary 3.3]. 

On a historical note, the road from the seminal article of Lee and Yang [16] on the Ising 
model (i.e. /xq = (<5i +#_i)/2 in (12. 4|) ) to the treatise |18| of Lieb and Sokal spanned almost three 
decades, with important contributions by Asano [2], Suzuki |29| and Griffiths [8 J (all concerning 
discrete spins), Ruelle's proof of a more general zero theorem |23| for Ising spins (generalizing 
a contraction method first introduced by Asano [3]), and the work of Simon and Griffiths |27| . 
which, among other things, establishes (|2.4p for d/j,o(o~) = exp[— acr 4 — ba 2 ]do~, with a > and 
6 6 I. (This particular measure indeed satisfies fj2 . 3[> . see |19|. Example 2.7]. It arises in the 
lattice approximation to the (1^)2.3 Euclidean field theory, cf. [25], Chapters VIII and IX, for 
which the Lee- Yang theorem is shown to hold [27} Theorem 6].) 

One-component models with more general interactions 

Next, we consider Hamiltonians of the form 

fljw = E - E h * a * = H 1'°M - E - *)> ( 2 - 6 ) 

XcA : \X\>2 x£-A x£-A 

for one-component spins o~ x G R, where we have added a constant linear in {h x } for later 
convenience. Theorem [2] implies that if the Boltzmann factor Bp(o~iC) = exp [ — /3H^ ,0 (o~a)], 

viewed as a function on C' A ' for fixed j3 > 0, belongs to V a +(^+ )> f° r some a > 0, then the 
partition function Zp^({h x }x&A) corresponding to the Hamiltonian (|2,6p satisfies 

Zp,A{{h x }xeA) + 0, for h x G H+, x G A, (2.7) 

for any Lee- Yang measure (j,q as defined above Theorem [2] 

If hq = (5 1 + <5_i)/2 (Ising spins), reasonably explicit results are known. The partition 
function is then seen to be a multi-affine polynomial in the activity variables z x = exp[— 2(3h x ], 
x G A: _ 

ZfiAM*eA) = E E *(P) II z *> ( 2 - 8 ) 

XcA x&X 

where Ex(/3) = exp [ — (3H^°(ax \a x = — l,x G X, = l,x G A\X)1 . The Lee- Yang property 
(12. 7|) then asserts that Zp^{{ z x}x&A) does not vanish whenever \z x \ < 1, for all x G A. If Ha'° 
is invariant under o"a; 1— > — x G A ("spin-flip" symmetry), the coefficients satisfy Ex ((3) = 
Ea\x(P)j f° r an A" C A. In particular, if z x = z, for all x, this yields Zp^(z) = z' A ' • Zp^z^ 1 ), 
for all z / 0, and therefore Zp t \(z) = implies that \z\ = 1. This is the Lee- Yang circle theorem. 
Ruelle's recent results \24\ Lemma 8 and Theorem 9] yield the following theorem. 
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Theorem 3. Let Z^\({z x } x ^\) be any multi-affine polynomial of the form (j2.8|) . with coefficients 
Ex{fi) = Ea\x(P) > 0; f or a M X C A, and all (3 > 0, that satisfies 

ZflAMveA) + 0, if \z x \ < 1 Vx G A, (2.9) 

for all (3 > 0. TTien ZgA * s ^ e partition function obtained by choosing <E>, in (|2,6p . to 6e an /smg 
ferromagnetic pair interaction, i.e., &({x,y})(cr x ,o~y) = —J X yO~ x o~y, for some J xy = J yx > 0, and 
&(X) = 0, otherwise. 

This is a converse to the Lee- Yang theorem for Ising spins. Ruelle also shows that, for 
general interactions <I> in (|2.6p with spin-flip symmetry, the Lee- Yang theorem can only hold at 
sufficiently low temperatures. (We refer the reader to [24] for a precise statement, and to |14j 
for a concrete example.) To further illustrate this point, we consider the Hamiltonian (|2.6p with 

#*'°(<7a) := — ^2 Juo~ U , where a u = o~ x and Jjj G E, 

UcA: \U\=2,4 x&J 

which has a pair and a four-spin interaction, and is invariant under a x h-> — (jj, x G A. The 
corresponding partition function Z^a^a) is given by (|2 . 8[) . Using the identities JT^g^l icr^) = 
Sxct/^-'-^)'^'' 7 ^' f° r an y ^ C A (the term corresponding to X = is understood to be 1), 
which imply 

PuM := + cr x ) + - o x ) = 2'~\ u \ £ 

X even 

one deduces that -ff * ° ma y be expressed in terms of pu{o~u) rather than a u . Indeed, 

■"A =~Jq- Ju-Pu{vu), where Ju = < ^ , , , , 

,7 CA ^| =2) 4 I 2J t/- 2 Ll/ C XcA:|X|=4^, lf|C/|-2 

and Jo is an irrelevant additive constant, which we may neglect. We note that, if a\ is the 
unique spin configuration such that {x G A : a x = —1} = X, for some given subset X of A, the 
quantity pu(o~u) vanishes unless U C X or U C A\X, in which case pu(o~u) = 1- It follows that 
Z/3 : a( z a) = ^2xcA E x (P)z x = (Z^*Z^ 4 ]j (z\) (* is called Schur-Hadamard product), 

where Z^ A (z A ) = ExcA E x\^ X ^ k = 2,A, with coefficients Ep \f3) = Uucx or ucA\x ^ h , 

\U\=k 

for all X C A. This *-product representation of Zp \ is very useful, because it essentially allows 

one to consider Z^ A , k = 2,4, separately, as we will now see. By [24, Example 7(d)], zf\ 
satisfies (|2.9p at temperature f3 > whenever 



PJu = %PJu > In 2 (or Ju = 0), for all U C A with \U\ = 4, (2.10) 
i.e., at sufficiently low temperature, for given Jjj > 0, |C/| = 4. Similarly, one obtains from |24[ 

J /3,A 



Example 7(a)] that zf\ fulfills {211]) if /3Jt/ > 0, for all U C A with \U\ = 2, or, equivalently, if 



Ju > J Xi for all J7 C A with \U\ = 2, (2.11) 

X: C/cXcA and |X|=4 

which implies that Jj; > 0, for \U\ = 2, and also requires the four-spin interaction to be suitably 

small as compared to the pair interaction. Finally, since Zp t A, = Za ^ * Z^\, Proposition 2(c) in 
|24| (see also |2H Corollary 2.15]; a crucial ingredient of the proof is a contraction method first 
introduced by Asano |3 : , Definition 3]) yields that Zp \, subject to the constraints (|2.10p and 
(|2.1ip . satisfies the Lee- Yang property (|2.9p . which is consistent with Theorem [3] 
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Multi-component spins 

For models of classical iV-component spins, N > 2, we denote the spin variable at site x by 
&X = {o~ x '■ i = 1, ■ ■ ■ , N), (we use Latin superscripts to indicate the components). For N = 1, 
the "Lee- Yang region" HL|_ is the region where the Laplace transform of the single-spin measure 
does not vanish, cf. condition fj2 . 3[> . An example of a suitable generalization to N > 2 is to 
consider measures \i G T N (see above Theorem [2]) obeying 

J e hal d{jL(a) ^ 0, for Re/i/0, (2.12) 

i.e., the magnetic field h is assumed to point in the 1-direction. If, in addition to satisfying 
(|2.12p . /x is assumed to be rotationally invariant, then it follows |18t Proposition 4.1] that its 
Laplace transform fi(h), h = (h 1 , . . . , h ), belongs to Vq+(Qn), where Qjy = U iljj and 



For the Hamiltonian 



N 

n± = {h = (h\...,h N ) eC N : ±Reh x >J2\ hi \}- ( 2 - 13 ) 

i=2 

N N 

Ha = - E ~ EE^> ( 2 - 14 ) 

{i,i/}cA i=l a;GA i=l 

the following result |18| Corollaries 4.4 and 5.5] (see also [1] for the plane rotator model) is 
known. 



Theorem 4. Let J xy G R, for all x, y G A and i = 1, . . . , N, with 

N 

•4>El J U V2; ^eA, (2.15) 

i=2 

and assume that //q G is a rotationally invariant measure on M. N satisfying condition (12.121) . 
Then the partition function Zp^\({h x } X £A) , (5 > 0, corresponding to the Hamiltonian (|2.14p does 
not vanish whenever h x G f2]y, /or a// x G A. 

The "Lee- Yang" region Qjj in Theorem |4] where the partition function Zp^[\h x } X £\\ is 
non-zero, can usually be further enlarged; see for example |18^ Proposition 4.1]. The domain 
([2.13|) considered here will suffice for our purposes. In particular, the possibility to include small 
but non-zero transverse fields h x , i = 2,...,N, x G A, will prove useful in what follows; see 
Section [3] 

For plane rotator models, i.e., N = 2 7 Theorem H] is "optimal" in the sense that J xy > \J xy \ is 
a sensible generalization of (|2.2p . For N > 3, the above result is not entirely satisfactory. Indeed, 
the constraint (I2.15P is a ferromagnetic condition that forces the interaction to be anisotropic 
and hence entails an explicit breaking of the 0(A^)-symmetry. Ideally, one would like to replace 
that condition by the more natural constraint 

Jl v > max Ijlj- (2.16) 

This is possible for the classical Heisenberg model (N = 3 and fio the uniform distribution on 
the unit sphere). See Proposition [6] below. 
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Quantum spins 

Quantum lattice systems for which the Lee- Yang theorem is known to hold include rather general 
(anisotropic) Heisenberg models with suitable ferromagnetic pair interactions. Although we 
are primarily concerned with models of classical spins, we quote some results concerning such 
models for the sake of completeness and because they imply the improved result for the classical 
Heisenberg model mentioned above. We fix a region A CC Z d and a spin s G N/2. By a x , 
i 6 A, i = 1,2,3, we denote the i-th. component of the quantum-mechanical spin operator 
at site x, in the spin-s representation of su(2). The "hat" distinguishes these operators from 
classical spins. They act on the space (C 2s + 1 )®I A I anc l satisfy the usual commutation relations 
[<T^,(jy] = i8 x ,y£jkl& x °f generators of su(2). We consider the Hamiltonian 

3 3 

h a ,s = - e - EE^> ( 2 - 17 ) 

{x,y}C\ i=l zGA i=l 

for arbitrary s G N/2, and define the partition function by 

Q/3,A,s({M*eA) = (2s + l)- |A| -ix[^ V {- KjS {{h x } xeK ))], /3>0. (2.18) 

The following theorem is essentially due to Asano [3] (see also [30]) and |4, Theorem 3], where 
complex transverse magnetic fields are considered. 

Theorem 5. Assume that the couplings J xy in (|2,17p satisfy the ferromagnetic conditions J xy > 
\J xy \ and J xy > \J xy \, for all x,y € A. Then, for arbitrary s € N/2 and (3 > 0, the partition 
function Q/3 t A,s ({^xIxga) in (|2.18p is non-zero whenever h x 6 Q,^ (see (|2.13p ). for all x G A. 

For results somewhat more general than Theorem [5] see the references given above. Our 
formulation is tailored to our purposes. In particular, the following result is a corollary of 
Theorem [5j 

Proposition 6. Let H\ be the Hamilton function defined in (|2,14p . with N = 3, and let /j,q 
be the uniform measure on the unit sphere S 2 . If the ferromagnetic conditions J xy > \J xy \ and 
Jxy > \ Jxy\> f or a M x iV ^ A, are satisfied, the partition function ^,A({/iijieA)j (3 > 0, does not 
vanish whenever h x G Q^, for all x G A. 

Proof. We define Q^\ SC } {{h x } x ^A} to be the partition function corresponding to the Hamiltonian 

H { j[ e s sc) defined as in (|2.17|) . but with rescaled coefficients, J xy i— )• J xy /s 2 and h x i— > h x /s. From 
a general result on classical limits of quantum spin systems due to Lieb |17| . it follows that 

lim Q^I A sc ] ({h x } xeA ) = Zp A ({h x } xeA ), pointwise, for real h % x , i = 1,2,3 and x G A. (2.19) 

Moreover, since \Q^\ SC ] ({h x }xeA)\ < exp (/3\\H A ^ ||) (where ||-|| denotes the operator norm) 
and \\H A e s sc ^ || is bounded uniformly in s, one sees that, for arbitrary compact subsets K C C 3 I A I, 

sup \Q { ^({h x } xeA )\<oo. (2.20) 

sGN/2; {h x } xeA eK ' 

Assuming that A = {x\, . . . , x n }, we write hi = h Xi and consider the hypothesis 

(H k ) : ]im^Q^^j(h k ; {hj}j^ k ) = Zp^{h k ; {hjjj^k), pointwise, and Zp^{h k ; {hj}^) / 0, 
whenever h k G fl^, for arbitrary (fixed) hj G Q^, j < k, and hj G I 3 , j > fc, 
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for k = 1, . . . , n. We begin by showing (Hi). To this end, we first define the sequence of functions 
(/ s 1 jl ) s with fl^hl) = Q ( ^s( h b hl,h\,{hj} m ), for arbitrary (fixed) hj, h\ G R and hj G R 3 , 
j > 1. By (|2.19p . / 8 j converges pointwise for real h\ as s — > oo, and by (|2.20p . it is uniformly 
bounded on compact subsets of C. It thus follows by Vitali's Theorem |31| Section 5.21] that 
lim^oo fji(h\) = Zp^(h\; hj,hf,{hj}j^i), uniformly on compact subsets of C. Moreover, 
since /jj does not vanish in H+ by Theorem \5\ Hurwitz' Theorem [T[ p. 178] implies that 
Z pA (hl; hj,h\, {hj} m ) ± for h\ G H+. 

Next, we consider the functions fgi(hj) = Ql e ^° s [hj; h\, hj, {hj}j^i) , for arbitrary (fixed) 
h\ G H+, hj G K and /ij 6 I 3 , j > 1. By what we have just shown, fgi(hj) converges as 
s — > oo for all hj G R and the limit is non-zero. Hence, Vitali's Theorem, together with the 
bounds (|2,20p . yields that lim s _ > . 00 f^i(hj) = Zp^(hj; h\, hj, {hj}j^i) , uniformly on compact 
subsets of C; and, using Hurwitz' Theorem together with Theorem [5l one shows that this limit 
does not vanish for hj G {z G C : \z\ < Re h\}. Repeating this argument for the sequence 
fl^hf) = Q^% C }(h\] h\,hl,{h,j} m ), with arbitrary h\ G W+, hj G {z G C : \z\ < Reh{} and 
hj G R 3 , j > 1, one obtains (Hi). 

In order to complete the proof of Proposition [6l it suffices to show that (H^i) implies (H k ), 
k = 2, . . . ,n. The proof is completely analogous to the one of (Hi) (one considers subsequently 
three sequences of functions (/j^) , for i = 1,2,3). The only difference is that one invokes 

(H^-i) instead of ([2.19p to argue that the functions f] k (h\) = Qp^gij^'i h\, hj, {hjjj^k) > (f° r 
fixed hj G j < k, hi, hi G R, and hj G R 3 , j > k), converge pointwise, for real h\, towards 
a non-zero limit, as s — > 00. □ 



3 Analyticity of correlations 

Our aim in this section is to establish analyticity of (connected) correlations as functions of an 
external magnetic field h in the Lee- Yang regions M± = {h G C : ±Re h > 0}, for the models 
introduced in the last section. For the sake of clarity, we first focus on models of Ising spins. 
Extensions to other models satisfying the Lee- Yang theorem are outlined at the end of this 
section. Thus, we consider the Ising Hamiltonian H\ defined in (|2.ip . but with a homogenous 
magnetic field h x = h for all x, and we write H^ = HaL_q- F° r the purposes of this section, we 
further impose periodic boundary conditions and require the pair interaction couplings to satisfy 

Jxy = Jyx = Jx-yfi > 0, for all x ^ y, and J 0x < oo. (3.1) 

x^O 

The a-priori measure /jlq on R is assumed to be any even measure with compact support satisfying 
condition (|2.3p (e.g., fio = + <^i)/2, the example of Ising spins, or /xo = the 

example of a continuous spin uniformly distributed on [—1, 1].) The partition function at inverse 
temperature f3 > is denoted by Zp ,\(h). Theorem Q] holds for this class of models. 

We are interested in the connected correlation functions (cr xi ; . . . ; c Xn )\ a h (also called Ursell 
functions, or cumulants in a probability theory context), where n G N and {x\, . . . ,x n } C A, 
which can be defined as [2b\ Section 11.12] 



(a xi ; . . . ; o- Xn )\^ h 



a £l 8 . 8£n log ^ exp E g ^]) w ) 



(3.2) 



where ( • )a,/3,/i denotes a (finite- volume) thermal average, i.e., an expectation with respect to 
the probability measure (Zp^(h)) e ~^ H ^^ Ylx&A c ^/ i o( cr x)- Our aim is to prove the following 
result. 
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Proposition 7. Under the above assumptions, for arbitrary sites x±, . . . , x n , n G N, and for all 
/3 > 0, i/ie thermodynamic limit 

((T xl ;...;a Xn )% h := lim (a xl ]...-,a Xn )%g h 

of the connected correlation (o~ xl ; . . . ; <J Xn ) A a ^ exists and is an analytic function of the magnetic 
field h, for Re h ^ 0. 

Proof. Since (3 > can be absorbed into a redefinition of J xy and /i, there is no loss of generality 
in setting (3 = 1. We consider the modified partition function 

/n 
exp{ -H A (a A ) + Y J [h + ^2e a e ika X y x }Yldfi (a x ), (3.3) 

with e = (ei,...,e n ) £ C n and fc a G A*, for all a = l,...,ra, where A* denotes the dual 
lattice of A. Assuming that Re(/i) > 0, setting h x := h + *YTa=\ e c* e%ka X anc ^ noting that 
Re(^x) > Re(/i) — Yla=i I £ q;|j Theorem [1] is seen to have the following corollary: 

n 

Z A (h,e)^0 if Re(/t) > J^|e«|. (3.4) 

a=l 

More generally, Z A (h,e) does not vanish whenever |Re(/i)| > Y^a=i \ £ oi\: because Z A ({h x } x ^ A ) = 
Z\({—h x } x£ \), by symmetry, since jiQ is even. In the following, we may therefore assume that 
Re(/i) > 0. Defining the convex domain 

n 

D = Uh,e) G C" +1 : Re(h) > ^ |e«|} C C n+1 , (3.5) 

a=l 

it is an easy exercise to check that Z A {h,e) has an analytic logarithm in D, so that 

/a(M) == log Z A (h,e) (3.6) 

is a well-defined, analytic function of (h, e) in D: (we choose a determination of log that is real, 
for h > and £ = 0). Henceforth, we fix some wave vectors k%, . . . , k n and consider a family 
of rectangular domains, A, with the property that A* 3 k a , V<x We then consider the analytic 
function 

TZ A (h,e) :=Z A (h,e)W := exp [/ A (M)L for (/i, e) G D. (3.7) 

We note that 7Z A is uniformly bounded in A on arbitrary compact subsets of D: abbreviating 
-HxeA [ h + Tl=i e o t e ika ' x ]o- x by H h ^ £ (<j A ), definition (J23]) implies that, for all (h,e) G D, 

\TZ A (h,s)\ = |Z A (/i,e)| 1/|A| < (He-^A+^A' 6 )^) 1 /^! < expflAI-^II^IU + H^IU)]. 
Clearly, 

n 

Halloo <c|A|-(h + ]Tm 

a=l 

for some constant c > 0; (here, we assume for simplicity that the support of the measure /iq 
is compact). Moreover, denoting by the ferromagnetic pair interaction, i.e., &({x,y}) = 
—Jxy°~xO~y, for all x ^ y, $>(X) = else, we have (see for example Section II. 3], and recall 
the definition of ||-|| below (|1.2p ) that 

II^aIIoo < c'|A| • p|||, 
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where c' > may depend on the dimension d of the lattice and takes into account the periodic 
boundary conditions imposed at dA. Notice that this last estimate holds for arbitrary translation- 
invariant interactions <£. The last two estimates imply that \TZ^(h,e))\ < c(h,e), where c(h,e) 
depends continuously on (h, e) and does not depend on A. Thus, for any compact K C D, 

sup \n A (h,e)\ < oo. (3.8) 

A;(h,e)eK 

With uniform bounds at hand, we may study the thermodynamic limit A /* Z d . With the 
help of a cluster expansion (see for example |32| ) at large magnetic fields, one shows the existence 
of the limit 

foo(h,e)= lim f\(h,e) (3.9) 

A/>Z d 

at any point in the (large-field) regime 

Re(/i) > h , \lm(h)\<6, \e a \ < 6, Va = l,...,n, (3.10) 

for sufficiently large ho and small 5 > 0; (had we kept the /3-dependence explicit, ho and 5 would 
depend on f3). Similar conclusions hold for 1Z\. Let S denote the subregion of D determined by 
the constraints (|3.10p . At any point (h,e) G S, the limit TZoo(h,e) = lim A ^ Z d lZ\(h, e) exists 
and is finite, and 'R. 00 {h,e) = exp [foo(h, e)l, by (|3,7p . In particular, TZ^ / on 5. Moreover, 
since the subregion S C D is a determining set for D, and since the functions TZ\ are uniformly 
bounded on compact subsets K of D, by (I3.8p . it follows from Vitali's Theorem |31|, Section 5.21] 
that TZ\ converges everywhere in D, as A /* Z , uniformly on any such K, and from Weierstrass' 
Theorem [1] p. 176] that the limit is analytic in D. 

Next, we show that TZoq vanishes nowhere on D. To this end, we first consider the functions 
lZ\{h,0) = exp[f\(h, 0)]; f\(h,0) is known to have a (pointwise) limit for real h |26| Section 
II. 3]. Hence TZoo(h,0) ^ 0, for arbitrary real h > 0. It follows from Hurwitz' Theorem [U 
p. 178] that TZ A (h,0) ^ 0, for all h € M+. Let (h,e) be any point in D. We show that 
TZ-oo(h, e) 7^ by repeated application of Hurwitz' Theorem: Define g\(z) = IZ^ih, z,0, . . . ,0), 
which is analytic and converges uniformly on compact subsets of D% := {z G C : \z\ < Re(/i)} 
towards g^z) := 7£oo(/i, z, 0, . . . , 0). But 5^(0) 7^ 0, hence by Hurwitz' theorem, g^z) vanishes 
nowhere in D\. In particular, ^(ei) = 7£oo(/i, e\, 0, . . . , 0) 7^ 0. Repeating this argument n times 
(in the A;-th step, the functions are g\(z) = lZ&(h, e±, . . . , £fc_i, z, 0, . . . , 0) and their domain of 
definition is = {z € C : \z\ < Re(/i) — Yli=i l e «l})> we obtain the desired result, namely that 

Having established that TZoo is nowhere vanishing on the convex domain D, it follows that 

foo(h,£l,. •-,£«) := log [Hoo(h,£i, . . .,£n)] 

is well-defined and analytic in D, and thus analytically continues /oo previously defined in (|3.9p 
on the subregion S of D determined by the constraints (|3.10p . Moreover, the functions f\ defined 
in (|3.6p converge towards /oo, uniformly on compact subsets of D. Similarly, any derivative of f\ 
with respect to variables (h,e) converges towards the corresponding derivative of /oo, uniformly 
on compact subsets of D. This follows from Cauchy's integral formula for polydiscs. 

Next, we consider correlation functions. Let the (discrete) Fourier transform of spins on Z rf 
and the corresponding reverse transformation be denoted by 

°k = Yl eik ' X(T xi keA*, and a x = — e~ ix ' k a k , x G A. 
xeA ' ' keA* 
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Observing that Z\(h,e) in (|3.3p can be rewritten as Z\(h,e) = Z\(h) ■ (^e^- ,a = 1 £a<Tfe a \ where 
( ' }a,/i = ( ' )h,p=l,hi it follows from definition (|3,2p that 



a" log [z A (M)] 



9ei • • • 9e„ 



=£n=0 



{a kl ;...;a kn ) c Kh 



(3.11) 



where the constraint Xla=i = [0] (with [0] = aL d for A = A(a) = 1, d /aL d ) follows from 
translation invariance, which holds because we have imposed periodic boundary conditions. With 
(|3.6p . we thus obtain from (|3. 1 1|) 



d n fA(h,e) 



de\ ■ ■ ■ de r 



(3.12) 



Proposition 7 then follows upon letting A /* using that the derivative on the left-hand side 
has a well-defined limit analytic in h in the region H+, and, subsequently, Fourier-transforming 
back to position space. (This does not affect analyticity in h, because all integrations in fe-space 
extend over a compact set.) □ 

Note that Proposition [7] continues to hold for unbounded single-spin measures fiQ satisfying 
suitable decay assumptions at infinity (this requires a somewhat more careful analysis). In 
particular, this is of interest in applications to field theory (see the end of Section 0] below). 

Next, we discuss generalizations of Proposition [7J to other models satisfying a Lee- Yang 
theorem. For one-component spins, we consider, for example, the Hamiltonian (|2.6|) . with a 
uniform external field h turned on, and = (<5-i + S\)/2. We assume that the interaction 
$ has spin-flip symmetry and that it is such that the Lee- Yang theorem holds at some inverse 
temperature /3 > 0. Note that this requires /3 to be sufficiently large, depending on $, if $ is not 
the ferromagnetic Ising interaction; see the discussion following Theorem [3] above and references 
therein, in particular J24J. Then {a xi ; . . . ; (T Xn )% h 1S analytic in h in the regions ELj.. The above 
proof is still applicable: the uniform bounds (|3.8p hold for general <I> with |||3>||| < oo, but the 
large-field cluster expansion, c.f. (|3,9p . must be modified slightly. 

The analyticity results of Proposition 7 can also be extended to certain A^-component models 
with N > 2. Consider the Hamiltonian (|2.14p and assume that condition (|2.15p holds, that the 
a- priori measure hq satisfies the assumptions of Theorem [Hand, in addition, that supp(/io) C 
is compact. Then the corresponding partition function satisfies 



N 



Zp,k{{h x } xeA ) ^ 0, if Re(hl) > ^ |, for all x G A. 



(3.13) 



i=2 



We assume that h x = (h, 0, . . . , 0), for all x in (|2.14p . To conclude analyticity of (c*^; . . . ; c^™ )o h 
in ELf., equation (|3.3p must be modified to read 



(h,e 1; 



exp 



P 



(7, 



n 
a=l 



E 

xGA 



c u x 



n d v°( c 

xeA 



.)■ 



Note that the conditions in (|3.13p are satisfied if Y2a=i l £ al ^ Pe(/i), which is the same constraint 
as for N = 1, c.f. ([3.4|) . The partition function Zgj^(h,e) is thus non-vanishing and possesses 
an analytic logarithm on the region D C C n+1 defined by ([33]). Hence, our proof of Proposition 
[7J applies in this case as well. In particular, the results hold when /j,q is the uniform measure on 
the unit sphere S N ~ l , and, for N = 3 (classical Heisenberg model), (|2.15p can be relaxed to the 
more natural ferromagnetic condition (|2.16p , by virtue of Proposition [6l 
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We conclude this section by discussing joint analyticity properties of correlations in h and 
/3, which follow from the proof of Proposition [7] above, using an idea of Lebowitz and Penrose; 
see |12^ Sections III and IV]. For simplicity, we restrict ourselves to the original framework of 
Proposition [7] but the following arguments apply to the other models discussed in the previous 
paragraph as well. We claim that (in the setting of Proposition [7J) , 

(a xi ; . . . ; o"x„}s h * s j° m tly analytic in f3 and h for (3 > and Re h > 0, (3-14) 

i.e., in some (complex) neighborhood of (0, oo) x H + in (J3, /i)-space. Indeed, this can be seen as 
follows. Let Zjy(P,h,e) be the modified partition function (|3.3|) . but with a factor of f3 inserted 
in the exponent, and the functions f\(f3,h,e) and TZ\(/3,h,e) be formally defined in terms of 
ZtSfit h, e) as in (|3.6p and (|3.7p . respectively. For arbitrary /3q > 0, by virtue of a suitable cluster 
expansion, one may select h sufficiently large and 5 = 5(/3q, h) > such that 

T^oo(l3,h,£) := lim TZ\(f3,h,e) is jointly analytic in the 

^ (315) 
variables (j3,h,e) on the polydisc Dp (5) x D^{8) x f_D (#)) , 

where D z (r) = {z' G C : \z' — z\ < r} denotes the open disk of radius r centered at z. For 
arbitrary ho G H+, let K C H+ be a closed rectangle containing both ho and D^(5) in its 
interior, and let / = Dp {8) Pi M+ be the (real) diameter of the disk Dg (8). The proof of 
Proposition [7J yields that 

7Zoo((3,h,£) is analytic in (h,e) on K x (Dq(5)) , for every (3 a I, (3.16) 

provided <5 > is sufficiently small. Moreover, it follows from (|3.8p that 



sup |^oo(/3,/i,e)| < oo. (3.17) 

Together with (|3.15|) . (|3.16p and (|3.17|) . the Malgrange-Zerner theorem [28 1 Theorem 2.2] (see also 
the lemma in Section III of [12]) then implies that 7£oo(A h, e) is jointly analytic for /3 G /, h G K, 
and e a G Dq(5), a = 1, ... ,n (i.e., in some neighborhood of this set in C 2+n ). In particular, 
TZoo((3,h,£) is jointly analytic in a polydisc around ((3o,ho,0), and we may further assume that 
7£oo is nowhere vanishing within this polydisc, by continuity (since 7£oo(A)j ^0) 0) 7^ 0, as shown 
in the proof of Proposition [7|). Thus, = logT^-oo is analytic on this polydisc, and so is 



E"=l k a = 



by ()3.12p . Since (3q > and ho G H+ were arbitrary, ()3.14p now follows as before upon Fourier- 
transforming back to position space, which does not affect analyticity in /3 and h. 



4 Applications 

In this section, we discuss various applications of Proposition [7] and of Eq. (|3.12p . 

Our first application concerns the independence of our analyticity results of the choice of 
boundary conditions; (recall that, in the above proof, we have imposed periodic boundary con- 
ditions). We wish to show that Proposition [7] continues to hold for all boundary conditions, b, 
for which the Lee- Yang theorem holds - this includes, in particular, free boundary conditions - 
and that the thermodynamic limit of correlations is unique for this class of boundary conditions, 
provided h belongs to the Lee- Yang region H_|_. 
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For the sake of simplicity, we sketch our arguments in the setting of Proposition [JJ but the 
following conclusions continue to hold for any of the Lee- Yang models (<3?,//o) mentioned above 
for which (an analogue of) Proposition [7J has been shown to hold. Thus, let Ha be the Ising 
Hamiltonian ([2.ip with h x = h, for all x, and pair couplings J xy satisfying (|3.1I) . and let /j,q be 
an even measure with compact support satisfying condition (|2.3p . Using a cluster expansion of 
the correlations at large magnetic fields, one proves (see for example |26| Theorem V.7.11]) that 
the model has a unique equilibrium state in a region fig C H_|_ defined by the constraints 

Be(h)>ho(p), \lm(h)\ < e(P), 

for some large ho > and some e > (both depending on f3). Denoting by ( • )as/i6 the 
(finite-volume) thermal average corresponding to a boundary condition b, it follows that 

A %]A| ( ^ i; --- ; ^ J W (41) 
exists and is independent of 6, for any h € fig. In particular, the limit agrees with 

}i™TT\(°ki-, ■■■;&>! Jl,B,h> ( 4 - 2 ) 

A/*Z d \A\ 

for h £ fig, where, in the latter correlations, periodic boundary conditions are imposed, as in 
Section 3. Note that the correlation in (|4.2p vanishes unless ^a=l = M ano - i s fi m t e if 
the latter condition holds. Assuming that the boundary condition b is such that the Lee- Yang 
theorem holds for h € H+, we may use Eqs. (|3.1ip and (I3.12p and then apply the Lee- Yang 
theorem to the partition function and the free energy of the model with boundary condition 
b imposed at dA to show that the correlations •••> &k n )x g hb are uniformly bounded 

and analytic in h on H + . Thus, multiplying ^{^ki', ■■■'■> ^k n )\^^ h by exp(-z^«=i k a ■ x a ) 
and integrating over the surface defined by the equation Y^a=i = [0], we find, using that the 
integration domain is compact, that 

lim (a Xl ]..r,a x Jlg hb =(a Xl -,...-,a Xn ) c g h , for all h G H + , (4.3) 

which proves analyticity of the thermodynamic limit of correlations with boundary condition b 
on the entire half-plane H + . 

Next, we consider the magnetization ((T x )p,h, which, by Proposition [7J is an analytic function of 
h in HL|_, for any f3 > 0. Since (3 is fixed in the sequel, we omit it from our notation. We propose 
to show that 

(cr x )h is a strictly positive, increasing, concave function of h > 0. (4-4) 

Note that this yields a well-known bound on a critical exponent for the magnetization as a 
function of h, (with j3 = /3 C , the critical inverse temperature). As a preliminary step, we show 
that (<To; Gx)h is decreasing for h > 0. We recall that 

d 

— (a Xl ; . . ■;<T Xn )% th = ^(er Zl ;. • ■ ; <r Xn \ °z)%, h i ( 4 - 5 ) 

Applying this identity for n = 2 and using the GHS-inequality (see [9] for Ising spins, |27| for 
more general /^o), we see that 

d 

z£A 
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for all h > 0. Letting A /* Z , one obtains that ^(o"o! < 0, for all h > 0. Indeed, this 
follows from Cauchy's integral formula, using that (oo; <r x )% h tends towards its infinite-volume 
limit uniformly on compact subsets of H_|_. (The latter claim follows from the proof of Proposition 
[TJ invoking Vitali's theorem, it suffices to establish uniform bounds, supa ; /J(0o> a x)\ /J < oo, for 
h belonging to an arbitrary compact subset K of H_|_, which, in turn, follow immediately from 
(HEED, iD, and the fact t hat (^;^>A,h = l A l _2 EfceA* e ik - x (a k ;a- k ) c Ath .) 

Returning to (|4.4|) . the identity (|4.5|) and the GHS-inequality imply that d 2 (a x )h/dh 2 < 0, 
for /i > 0. Similarly, monotonicity of (a x )h, for /i > 0, follows from (|4.5p and the FKG-inequality 
(which holds since the interactions are ferromagnetic). It remains to show that {cr x )h is positive 
for h > 0. Indeed, {a x )h > 0, for all h > 0, follows from Griffiths' inequality. If {cr x )h vanished, 
for some h > 0, then {<J x )y = 0) f° r an < /t' < by monotonicity, which is impossible, because 
the zeroes of ((T x )h form a discrete subset of H + . Similarly, the first derivative of {cr x )h in h is 
strictly positive and the second derivative has, at most, a discrete set of zeros, for h > 0. (Note 
that, since {<J x )-h = — {&x)hi one derives corresponding properties of the function (cr x )h, h < 0, 
from flOJ) .) 

Next, we consider the mass gap, m((3,h) (inverse correlation length), defined as 

m (^ h ) = era u\ = - limsu P T-T lo g IWo'^x)^,^- (4.6) 
slPi^-J |aj|->oo \ x \ 

We restrict our attention to the generalized Ising model considered in Proposition [7] with finite- 
range interactions (i.e., J xy = whenever \x — y\ > R, for some R > 1), but the following 
discussion applies to more general spin models obeying a Lee- Yang Theorem and to Euclidean 
A(/> 4 -field theories with non-zero "external field" h, in two and three space-time dimensions; see 
|12| . |13| . and |10| . The mass gap satisfies 

m(/3, h) > 0, for all h G H + and /3 > 0, (4.7) 

i.e., the two-point function (o"o; & x )p ^ exhibits exponential clustering. The inequality (|4.7p can 
be shown as follows. By Proposition [71 (<Jo', a 'x)h is analytic in H + , hence log | (ctq; CTx)^ is a 
subharmonic function of h € H+, see |10^ Theorem A. 3]. Using a cluster expansion, one shows 
that m is positive on an open subset of H + corresponding to sufficiently large Re h; (actually, 
positivity on some smooth arc in M + would suffice). It then follows from subharmonicity (see 
[13, Lemmas 1 and 5]) that m is positive everywhere on H_|_. 
Next, we remark that 

m(/3, h) is an increasing function of h > 0, for all f3 > 0. (4-8) 

This follows from definition (|4.6p . using monotonicity of the logarithm. For, we have already 
shown that the two-point function (<Jo]<J x ) C p h is decreasing in h, for h > 0. We also wish to recall 
a bound on the critical exponent 5 describing the divergence of the correlation length £ = m _1 , 
as h \, 0, at the critical inverse temperature /3 C ; (see |10| for a more general result that extends 
to Euclidean A0 4 -field theory): 

£(/3 c , h) ~ h- s , with 6 < 1. (4.9) 

This follows by showing that m(/3 c ,h) > c ■ h, for some positive constant c and all sufficiently 
small h > 0. The latter is a consequence of Theorem A. 6 in |10] . using the fact that the functions 
— log | ((To! a x)p c h\ are superharmonic in H + and that, given any (real) ho > 0, these functions 
are bounded away from by a positive constant (uniform in h and x), for all h > ho and all 
sufficiently large x, which follows from (|4.7p and (|4.8p . 
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Finally, we mention a generalization (alluded to, above) of PropositionUJto some TV-component 
(N = 1,2,3) Euclidean A|0|^-field theories (<f> = ((f) 1 ,- ■■ ,<ft N )) in d = 2,3 space-time dimen- 
sions, with periodic boundary conditions. The correlation functions of lattice spin systems are 
replaced by the Schwinger (Euclidean Green) functions, 

S L,a,h(%l j ) • • • j %m Oin) , 

of the properly renormalized lattice field theory on a lattice n Ax,, with periodic boundary 
conditions imposed at dAi, where d = 2,3, a > is the lattice spacing, Al = [-L/2, L/2] d is 
a cube in Z„ with sides containing ^ sites, and the arguments Xi,ai, i = 1, . . . ,n, stand for the 
field components (j) ai (xi). Formally, the Schwinger functions are given by 

J> Ql(a:i) • • • (f> a -(x n )e- A ^M^D<j> Lta 



J e -A L , a , h (4>)Dcj) L a 

where A^^^ is the Euclidean action of the theory, which is identical to the Hamilton function 
with periodic boundary conditions of the corresponding classical lattice spin system (with nearest- 
neighbor couplings, and in an external magnetic field h), but with coupling constants that depend 
on the lattice spacing a in such a way that the continuum limit a — > exists; see, e.g., [25], |20| . 

Combining results in [25j and |28j . for d = 2, and in [20j and [5], [28], for d = 3, one can 
prove an analogue of Proposition [JJ for h € H+ . This is accomplished by first proving existence 
and Euclidean invariance of the limits 



lim ]3mS Lah (xi,ax,...,x n ,an), (4.10) 

for Re/i large enough and Imh small enough. The Lee- Yang theorem (see [271 Theorem 6]) and 
uniform bounds on the analogue of the free energy, see (|3.6p . discussed in [25j (d = 2) and in [7] 
(d = 3) then imply analogues of Eqs. (|3.8p and (|3,12p that can be used to prove bounds on the 
Schwinger functions, integrated against test functions on momentum space, that are uniform in 
a and L and yield analyticity in h on M + . As a consequence, the limiting Schwinger functions in 
(|4.10p exist, are Euclidean invariant and analytic in h, for h G H+. A more detailed discussion 
of these arguments goes beyond the scope of this note. 

The results mentioned here are of interest in an analysis of phase transitions accompanied 
by spontaneous symmetry breaking in A|0| 4 — theory in d = 3 space-time dimensions; see |6, 
Section 4]. 



Acknowledgements. We thank David Ruelle for informing us about ref. |14) prior to publication 
and for several helpful discussions. We are grateful to Barry Simon for helping us to trace some 
references. The senior author thanks Elliott Lieb for all he has taught him and for his friendship. 



References 

[1] Ahlfors, L. V.: Complex Analysis, 3rd ed., McGraw-Hill (1979) 

[2] Asano, T.: Generalized Lee-Yang Theorem. J. Phys. Soc. Jap. 25, pp. 1220-1224 (1968) 
[3] Asano, T.: Theorems on the Partition Functions of the Heisenberg Ferromagnets. J. Phys. Soc. Jap. 29, 
pp. 350-359 (1970) 

[4] Dunlop, F.: Analyticity of the Pressure for Heisenberg and Plane Rotator Models. Commun. Math. Phys. 
69, pp. 81-88 (1979) 

[5] Feldman, J. and Osterwalder, K.: The Wightman Axioms and the Mass Gap for Weakly Coupled (0 4 )3 
Quantum Field Theories. Ann. Phys. 97, pp. 80-135 (1976) 



15 



[6] Frohlich, J., Simon, B., and Spencer, T.: Infrared Bounds, Phase Transitions and Continuous Symmetry 

Breaking. Commun. Math. Phys. 50, pp. 79-95 (1976) 
[7] Glimm, J. and JafTe, A.: Positivity of the (4> 4 )3 Hamiltonian. Fortschritte der Physik 21, pp. 327-376 

(1973) 

[8] Griffiths, R. B.: Rigorous Results for Ising Ferromagnets of Arbitrary Spin. J. Math. Phys. 10, pp. 1559- 
1565 (1968) 

[9] Griffiths, R. B., Hurst, C. A., and Sherman, S.: Concavity of Magnetization of an Ising Ferromagnet in a 
Positive External Field. J. Math. Phys. 11, pp. 790-795 (1970) 

[10] Guerra, F., Rosen, L., and Simon, B.: Correlation Inequalities and the Mass Gap in P(<j>)2- Commun. 
Math. Phys. 41, pp. 19-32 (1975) 

[11] Lebowitz, J.: Bounds on the Correlations and Analyticity Properties of Ferromagnetic Ising Spin Systems. 
Commun. Math. Phys. 28, pp. 313-321 (1972) 

[12] Lebowitz, J. and Penrose, O.: Analytic and Clustering Properties of Thermodynamic Functions and Dis- 
tribution Functions for Classical Lattice and Continuum Systems. Commun. Math. Phys. 11, pp. 99-124 
(1968) 

[13] Lebowitz, J. and Penrose, O.: On the Exponential Decay of Correlations. Commun. Math. Phys. 39, pp. 
165-184 (1974) 

[14] Lebowitz, J. L. and Ruelle, D.: Phase Transitions with Four-Spin Interactions. Available at arXiv: 1006.0959 
[15] Lee, T. D. and Yang, C. N.: Statistical Theory of Equations of State and Phase Transitions. I. Theory of 

Condensation. Phys. Rev. 87, pp. 404-409 (1952) 
[16] Lee, T. D. and Yang, C. N.: Statistical Theory of Equations of State and Phase Transitions. II. Lattice 

Gas and Ising Model. Phys. Rev. 87, pp. 410-419 (1952) 
[17] Lieb, E. H.: The Classical Limit of Quantum Spin Systems. Commun. Math. Phys. 31, pp. 327-340 (1973) 
[18] Lieb, E. H. and Sokal, A. D.: A General Lee-Yang Theorem for One-Component and Multicomponent 

Ferromagnets. Commun. Math. Phys. 80, pp. 153-179 (1981) 
[19] Newman, C. M.: Zeros of the Partition Function for Generalized Ising Systems. Comm. Pure Appl. Math. 

27, pp. 143-159 (1974) 

[20] Park, Y. M.: Lattice Approximation of the (A</> 4 — fj,(j>)3 Field Theory in a Finite Volume. J. Math. Phys. 
16, pp. 1065-1075 (1975) 

[21] Rodriguez, P.-F.: Around the Lee-Yang Theorem and its applications. Master's thesis, Department of 

Physics, ETH Zurich, supervised by Prof. Dr. J. Frohlich (2011) 
[22] Ruelle, D.: Statistical Mechanics - Rigorous Results, Benjamin (1969) 

[23] Ruelle, D.: Extension of the Lee-Yang Circle Theorem. Phys. Rev. Let. 26, pp. 303-304 (1971) 
[24] Ruelle, D.: Characterization of Lee- Yang Polynomials. Ann. Math. 171, pp. 589-603 (2010) 
[25] Simon, B.: The P(4>)2 Euclidean Quantum Field Theory. Princeton University Press (1974) 
[26] Simon, B.: The Statistical Mechanics of Lattice Gases, Vol. 1, Princeton University Press (1993) 
[27] Simon, B. and Griffiths, R. B.: The (0 4 ) 2 Field Theory as a Classical Ising Model. Commun. Math. Phys. 
33, pp. 145-164 (1973) 

[28] Spencer, T.: The Mass Gap for the P(</?)2 Quantum Field Model with a Strong external Field. Commun. 

Math. Phys. 39, pp. 63-76 (1974) 
[29] Suzuki, M.: Theorems on the Ising Model with General Spin and Phase Transition. J. Math. Phys. 9, pp. 

2064-2068 (1968) 

[30] Suzuki, M. and Fisher, M. E.: Zeros of the Partition Function for the Heisenberg, Ferroelectric, and General 

Ising Models. J. Math. Phys. 12, pp. 235-246 (1970) 
[31] Titchmarsh, E.C. : The Theory of Functions. 2nd ed., Oxford University Press (1939) 
[32] Ueltschi, D.: Cluster Expansions and Correlation Functions. Mosc. Math. J. 4(2), pp. 511-522 (2004) 



16 



